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For a doped antiferromagnet with short-range spin stripe correlations and long-range charge stripe 
order we find that the manifestation of charge order changes abruptly as a function of momentum 
along the Fermi surface. The disorder averaged local density of states is almost perfectly homoge- 
neous when integrated only over states which contribute to the "nodal" spectral weight whereas it 
displays long range charge stripe order when integrated only over states which contribute to the 
"antinodal" spectral weight. An effectively two dimensional isotropic nodal liquid can thus coexist 
with static charge stripes provided there is no static spin order. We also study commensurate spin 
and charge stripe ordered systems where the Fermi surface consists of a nodal hole pocket and an 
open "stripe band" section. Due to the stripe order the relation between hole density and size of a 
pocket will be reduced compared to a paramagnet by a factor of two for even charge period and four 
for odd charge period and we find an estimated upper limit on the area fraction of a hole pocket 
of 1.6% for charge period four and 4% for charge period five. We also discuss why electron pockets 
are not expected for a stripe ordered system and show that the open Fermi surface section may be 
electron like with a negative Hall coefficient. 



I. INTRODUCTION 

Lately it has become increasingly clear that an intrinsic 
property of the cuprate high-temperature superconduc- 
tors is the presence of nanoscale inhomogeneity of the 
low-energy electronic spectral weight in the doped Cu02 
planes. In particular this inhomogeneity seems to come 
in the form of stripes which are unidirectional spin and 
charge density modulations. From neutron scattering ex- 
periments for example there is evidence for a universal 
spin response in the hole doped materials that is con- 
sistent with dynamic spin stripe correlations.'^ The ten- 
dency for stripe formation in a doped antiferromagnet is 
well established and can be understood naively as an ef- 
fect of the competition between strong antiferrom agnetic 
correlations and delocalization of the doped holes .'^l^liHl 

Concurrently with the increasing evidence for stripes 
there is also a picture emerging in which there are 
electronic quasiparticles in the normal state and Bo- 
goliubov quasiparticles in the superconducting state in 
line with that expected for a Fermi liquid or a BCS 
superconductor.'^'-^ In particular, for underdoped mate- 
rials these quasiparticles appear to be confined in mo- 
mentum space to the so called nodal region of the Bril- 
louin zone, near momentum (7r/2, 7r/2), where the d-wave 
gap has a node in the superconducting state. The recent 
observations by scanning tunneling spectroscopy (STS) 
of a "stripe glass" in which there are static unidirec- 
tional electronic domains makes the problem even more 
intriguing.'^ In fact, this suggests that the two phenom- 
ena of stripes and nodal quasiparticles can coexist. A 
most basic question is how well defined two dimensional 
quasiparticles can exist in a system where there is one 
dimensional charge order. 

In this paper we consider a mean-field model of a sys- 
tem with long-range charge-stripe order and short-range 



spin-stripe correlations that we suggest may hold some 
clues to this puzzle. From Landau theory describing the 
low-energy long-distance physics it is known that spin 
stripe order alw ays implies charge stripe order but not 
the convers^^MI thus providing the logical background 
for the study of such a system. Assuming translational 
symmetry along the stripe direction allows us to define a 
momentum resolved local density of states (LDOS) which 
is integrated only over states with a particular longitu- 
dinal momentum. Although this may not be an experi- 
mentally measurable quantity it will give some essential 
insights into the significance of the lack of long-range spin 
order. 

What we find is an abrupt change in the manifestation 
of the charge stripe order at low energies as a function of 
the momentum along the Fermi surface. Quite strikingly, 
as shown in Figure [S] the spin disorder averaged LDOS 
at the Fermi energy when integrated only over the states 
which contribute to the nodal region spectral weight show 
very little evidence of the charge order and is nearly ho- 
mogeneous and isotropic, while the full LDOS including 
also antinodal, near (7r,0), states clearly display charge 
order. This should be contrasted with results for a sys- 
tem with both charge and spin order, as shown in Figure 
|3] for which the LDOS has a clear signature of charge 
order for antinodal as well as nodal states. The results 
thus suggest that an effectively two dimensional, isotropic 
nodal state can indeed coexist with static charge stripe 
order provided that there is no static spin order. We have 
considered long range charge order but the results are ex- 
pected to apply equally well for a system with only local 
stripe correlations provided that the charge order is effec- 
tively static and the spin fluctuations are fast compared 
to the experimental probe. 

Another basic question concerns the nature of the 
Fermi surface of the underdoped cuprates. Angle re- 
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solved photoemission (ARPES) experiments observe the 
famous nodal region "Fermi arc" which is most read- 
ily consistent with a large Fermi surface closed around 
the (tt, tt) point but where the spectral weight is sup- 
presse d by the so called "pseudogap" in the antinodal 
region PliE2ll3llil large Fermi surface picture has 
been challenged by the recent observations of high fre- 
quency quantum oscillations of the Hall resistivity in 
strong magnetic fields in YBa2 Cu30„ (YBCO) that seem 
to imply small Fermi pockets 1^^'^^'^ We show that ex- 
actly this distinction is found between a system with 
long-range spin stripe order and a system with only short- 
range spin stripe correlations. In simulations of commen- 
surate stripe ordered systems (Fig. [l]) we find a nodal 
region hole pocket together with a quasi-one dimensional 
Fermi surface section that are replaced by a large Fermi 
surface in a disordered stripe system (Fig. [5]) . It is known 
that a c-axis magnetic field can enhance stripe ordei'^^, 
thus providing a possible connection to the observations 
of quantum oscillations. If the high-field measurements 
indeed probe a stripe ordered state which the ARPES 
measurements generally do not our results may thus pro- 
vide an explanation for the apparent discrepancy between 
the two probes where the Fermi arc could be regarded 
as the remnant in the disordered state of a nodal hole 
pocket. 

In light of this we also consider the implications of the 
enlarged stripe unit cell and find that the hole density of 
a hole pocket will be decreased by a factor of two for even 
charge period and four for odd charge period compared 
to an estimate that assumes no long-range order. Based 
only on the stripe periodicity we can also set an estimated 
upper limit to the size of a nodal hole pocket of 1.6% of 
the full Brillouin zone for a period four stripe and 4% 
for a period five stripe with larger pockets merging into 
open sections. From the canonical relation between dop- 
ing and stripe periodicit}^ we may expect period five for 
"1/10" doping and period four for "1/8" doping. Quan- 
tum oscillations have been observed for doping close 1/10 
whereas for 1/ 8 doping the high field limit remains to be 
explored.'i^^For the period four stripes we expect either 
no oscillations, corresponding to only open orbits, or pos- 
sibly a smaller frequency of up to 450T, corresponding to 
the maximum pocket size. If a sharp distinction is found 
between 1/8 and 1/10 it would be a dramatic confirma- 
tion of stripe order. 

Although hole pockets may explain the frequency 
of quantum oscillations they cannot explain why 
the Hall coefficient (Rh) may be negative at low 
temperature.^^ This has lead to suggestions for the 
formation of ele ctron pockets due to broken translational 
symmetry ]^ ^ l ^^l However, the location of such electron 
pockets are in regions of the Brillouin zone (along (0, 0) 
to (tTjO)) where there is no evidence from photoemission 
of any substantial spectral weight and they must be con- 
sidered highly speculative. For a stripe ordered system, 
in particular, we argue that electron pockets arise only 
as an artifact of a mean-field type description that effec- 



tively ignores the interactions on a stripe. 

Here, instead, we show that the stripe band with open 
orbits may be electron-like thus providing an alterna- 
tive explanation for a negative Hall coefficient. The 
sign of Hh for the stripe band depends sensitively on 
the changes of the Fermi velocity over the Fermi surface 
which depends in turn on the band structure parameters 
as well as the character (strength and periodicity) of the 
stripe or der. The ap parent nonuniversality of the Hall 
coefficienillSIlSMESEil different materials is not unex- 
pected if stripe order plays a role. 



II. THE MODEL 

The model we consider is a tight-binding model on 
a square lattice in a static potential which couples to 
the local spin density and which may or may not have 
long-range stripe order. We can think of the poten- 
tial as a strong inhomogeneous magnetic field that is 
generated self-consistently from an interacting model 
such as the large-U Hubbard model in the Hartree-Fock 
approximation.!^ The Hamiltonian reads 

H = - tJ2 (4,.c.v + H.C.) 

{rr')(T 

- t' J2 (4,.c,.'. + H.C.) 

+ '^{~irv{x)ciy^,c,,y,. , (1) 

x,y,cr 

where Cr,a is the electron destruction operator at site 
r = {x,y) and with spin a — ±. We will also refer to 
the stripe perpendicular (x) and stripe parallel (y) direc- 
tions as x± (fcj_) and x\\ (fcy) respectively. The hopping 
is given in a standard fashion where (rr') indicates near- 
est neighbors, (rr')' next-nearest neighbors. We will use 
energy units such that t — 1 and we take t' — —0.3. For 
the Hubbard model with on-site interaction U nr^riri 
we would identify the quantity cr( — l)^y(x)/[/ with a 
self-consistent spin-density wave. (In a stripe ordered 
state there would also be a charge density wave gen- 
erated potential of strength USn/2, with 6n the local 
density variation, which will in general be quite negligi- 
ble compared to V.) In this paper we will not attempt 
to solve any microscopic model self-consistently, instead 
the Hamiltonian provides a simple means of simulating a 
system with stripe correlations by choosing appropriate 
potentials V{x). This type of model has been considered 
before in several publication^^fiisa however not to con- 
sider systems with only charge stripe order and neither 
has the local density of states been explored in any de- 
tail. We will always consider systems which are less than 
half-filled thus providing a most simple caricature of a 
hole-doped cuprate with stripe order or stripe correla- 
tions. We also consider only commensurate stripe order 
in which the stripe period is some integer multiple of 
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FIG. 1: (Color online) Low energy spectral weight of period 
four (a,b) and period five (c,d) bond-centered stripe ordered 
systems for stripe potential — 0.5 and doping given 

below. The left column shows the spectral weight as an in- 
tensity plot where the dashed lines are the Bragg planes of 
the stripe order at (37r/4,7r) and (47r/5,7r) respectively. The 
right column shows the full Fermi surface together with the 
centermost Bragg planes (dashed) and the first Brillouin zone 
(highlighted). The size of a pocket in b is 0.7% of one quad- 
rant (0 < fc|| < TT, < fcx < tt) and the actual hole density 
of the nodal pockets is also 0.7%, with a total hole doping 
including the antinodal stripe states of 17.5% {Ef — —1.15). 
The size of a pocket in d is 1.2% of one quadrant while the 
hole density of the nodal pockets is only 0.6%, with a total 
doping of 12.5% {Ef = —1.0). The boxes show the estimated 
upper limit to the size of a pocket as discussed in the text. 



the lattice spacing. Although fluctuating incommensu- 
rate stripe correlations are found experimentally, stripe 
ordered states always appear to be commensurate.l^ 



III. ORDERED STRIPES 

Let us first consider the case of a completely reg- 
ular periodic potential. To exemplify the results we 
consider bond-centered stripes with charge periodicity 
four and five. For the period four array we choose 
a potential Vi{x) — 0.5V4,o(a;) with 14,0(2; mod 8) = 
1,1,-1,1,-1,-1,1,-1 and for the period five array 
V5(x) = 0.5V5,o(a;) with V:^,o(a;mod 5) = 1,1,-1,1,-1. 
Note the size of the magnetic unit cell transverse to the 
stripes which is twice the charge period only for even pe- 
riodicities. Considering now also the antiferromagnetic 
(AF) (—1)^ modulation (Eq. [T]) along the stripes we 
have primitive vectors of the magnetic unit cell (7r/4, ±7r) 
for period four and (27r/5,0) and (0, tt) for period five. 
The modulation of the potential is basically the sim- 



plest possible consistent with bond-centered stripes. A 
self-consistent potential is expected to have a somewhat 
smaller amplitude on a stripe than between stripes be- 
cause of the higher hole density but such a consideration 
would not change the qualitative features of our results. 

The standard mean field treatment of an antiferromag- 
net gives a lower and upper Hubbard band separated by 
an indirect gap (for t' ^ 0) while the stripes in the form 
of TT phase shifts of the potential transfer some of that 
weight to midgap states. With moderate doping, the 
chemical potential will lie in this midgap stripe band. 
Depending on the magnitude of the potential and the 
doping there may be additional states at the Fermi en- 
ergy coming from the stripe modified lower and upper 
Hubbard bands. If the potential is large the lower band 
will be gapped away from the Fermi energy leaving only 
the midgap stripe states which as we will find have weight 
primarily in the antinodal region. If the potential is 
small the upper band will also cross the chemical po- 
tential, giving rise to electron pockets. For hole doping 
and considering as a guide the large-U Hubbard model 
such electron pockets are most likely unphysical and we 
make sure to choose to potential to avoid such states. 
Spectral weight in the nodal region most naturally arises 
from the lower band but does require a somewhat judi- 
cial choice of the magnitude of the potential. Whether 
the stripes are bond-centered or site-centered is however 
not important for the qualitative features of the Fermi 
surface. 

We diagonalize the Hamiltonian to get a spectrum 
uja with single-particle eigenstates ^'a(fc) (a includes 
band index and spin) that only have spectral weight 
at momenta connected by multiples of the recipro- 
cal lattice vectors of the stripe unit cell. The zero- 
temperature single-particle spectral function is simply 
given by A(fc,w) = J2c I*a(fc)pi5(w - Wa)- Figure [l] 
shows the Fermi surface (integrated over a narrow win- 
dow Soj ~ 0.02) for the two realizations of stripe poten- 
tials. For this figure we have chosen the Fermi energy and 
corresponding hole doping to give a simple Fermi surface 
with only two bands contributing, one band giving the 
nodal Fermi pocket and one band giving the quasi-one 
dimensional antinodal weight. The antinodal spectral 
weight is qualitatively not very sensitive to the Fermi 
energy, the Fermi surface only moves along the stripe di- 
rection fc|j as expected for a quasi-one dimensional band. 
The nodal weight on the other hand is very sensitive to 
the details of the potential and the doping, in particular 
for higher doping several bands may contribute and give 
superimposed nodal hole pockets with different sizes. 



Fermi surface 

We draw the following general conclusions from the 
calculations on ordered stripe arrays: (1) The size of a 
nodal hole pocket does not give the corresponding hole 
density without an a priori knowledge of the stripe order. 
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Naively, for a paramagnet with one pocket per quadrant 
of the full 1st Brillouin zone, one might infer a hole den- 
sity which is Ufi = 8A where A is the area fraction of a 
pocket. Here we find that the actual hole density is one 
half of this value Uh = A A for period four stripes and only 
one quarter rih = 2A for period five stripes, i.e. there are 
only two or one pocket respectively in the reduced first 
Brillouin zone. We have checked that these results gener- 
alize to longer stripe periodicities with the same distinc- 
tion between even and odd periodicity. Thus, only for 
even periodicity is the relation between hole density and 
the size of a pocket the same as for an antiferromagnet.'^ 

(2) The bulk of the hole density for this system is actu- 
ally in the quasi-one dimension al stri pe band which may 
dominate the charge transport .1221221 However, since the 
corresponding Fermi surface section does not have any 
closed orbit this band would not give rise to quantum 
oscillations. 

(3) The size of the Brillouin zone limits the size of a 
nodal pocket, larger pockets will merge to form extended 
Fermi surface sections (open orbits). The size limit de- 
pends on the shape and orientation of the pocket, how- 
ever we can make a general and conservative upper esti- 
mate by assuming that the extension of the pocket is of 
the same magnitude along both directions as indicated 
by the boxes in Figure [l] (b) and (d) . With this we es- 
timate the maximum area fraction (of the full extended 
BZ) to Ai,rnax = 1/64 ft! 1.6% for charge period four and 
A^^max = 1/25 = 4% for period five. 

The above estimate does appear to rule out a period 
four stripe ordered state in the system (YBa2Cu306.5) 
with hole doping p = 0.10 in which quantum oscillations 
were recently observed with a pocket of size 1.9%.^^ The 
size is however well within the limit for a period five 
state, which may be reasonable given the observations by 
inelastic neutron scattering of magnetic corre lations that 
are consistent with this stripe periodicitjHini Por charge- 
period five the size of the pocket would correspond to 
hole density of only 3.8%, while the remaining 6.2% hole 
density would according to this model reside in the quasi- 
one dimensional stripe band. 

The cyclotron mass of a pocket is given h^^^ ^ = 
— ^ « 2.61 %i where we have used a lattice spacing 

of 3.8A, t = OAeV, A is the area fraction of one quad- 
rant of the Brillouin zone and e is energy in units of t. 
A problem with these pockets in relation to the experi- 
ments is that the dispersion is quite steep which gives a 
cyclotron mass of around 0.3 electron masses for both the 
period four and period five stripes which is much smaller 
than the 1.9 reported experimentally. However, one may 
expect that strong correlations can significantly narrow 
the bands and increase the cyclotron mass. 



Electron pockets and open orbits. 

The Hall resistivity is negative under the conditions in 
which the quantum oscillations have been detected which 
from semiclassical transport theory implies electron like 
charge carriers instead of the expected hole like. This 
led to a recent suggestion that electron pockets, not hole 
pockets, in a stripe ordered system may be responsible 
for the oscillations .^23 In contrast to the hole pockets dis- 
cussed in this paper such electron pockets are centered 
on the stripes and may appear with a suitable choice of 
the potential. However, we would argue that they are 
unphysical as they correspond effectively to occupying 
states in the upper Hubbard band of the stripe, with the 
relatively low energy being an artifact of the mean field 
like treatment of the stripe. 

We can see this by considering a single site-centered 
stripe in the large U limit. Self-consistently this corre- 
sponds to a large staggered potential outside of the stripe 
center row giving rise to a decoupled tight binding chain 
on the stripe.'221 The problem with this is particularly 
obvious at half-filling of the stripe where the procedure 
gives a gapless electronic band instead of the appropriate 
insulator with a gap proportional to U. This discrepancy 
is related to the fact that the mean field treatment cannot 
accurately describe the effectively paramagnetic state on 
the site-centered stripe. Clearly the naive stripe potential 
misses essential aspects of the interactions on a stripe, a 
deficiency that will remain also away from half-filling and 
for a smaller stripe potential. 

The resolution to the negative Hall resistivity may lie 
inste ad in the charge transport related to the open stripe 
bands^^'^^'with the oscillations being due to hole pockets. 
For the transport properties of the quasi-one dimensional 
stripes it may be important to consider correlations^ as 
evidenced by the antinodal pseudogap but here we will 
do the standard semiclassical calculation. Following the 
presentation of On^2l!^ tJ^e sign of the Hall conductivity 
is given by the integral over the Fermi surface z ■ J dl x I 
where I = v^r-j: is the scattering length vector with Fermi 
velocity Vj^ and relaxation time r^. Figure [2^ shows con- 
stant energy curves of the "stripe band" of the charge 
period four system considered earlier (same parameters) . 
With increasing energy these evolve from electron pock- 
ets into open sections and into hole pockets with the Hall 
coefficient changing continuously from negative (electron- 
like) to positive (hole-like). (We have assumed a k inde- 
pendent relaxation time.) 

Figure |2]3 shows the area traced out by Fermi veloc- 
ities for the particular section marked by arrows as the 
Fermi surface is followed across the Brillouin zone where 
the velocities and v± corresponding to the longitudi- 
nal and transverse components respectively. We see that 
the outer part of the Fermi surface {k± ft 7r/4) gives a 
hole-like contribution whereas the inner part {k± ft 0) is 
electron-like. The particular Fermi surface displayed in 
Figure [l] has a hole- like open section but is close to the 
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FIG. 2: (Color online) Fermi surface sections (a) in the upper 
half of the first BZ of the stripe band for the charge period 
four stripe (same as in Fig. [l] a,b) for varying chemical po- 
tential. The three lower sections are electron-like and the rest 
are hole-like. As explained in the text the Hall coefficient is 
proportional to the area of the integral over the Fermi veloc- 
ity curve as shown in (b) for the section in (a) marked by 
arrows. There is a contribution from both a hole-like and 
an electron-like loop. The magnitude of the longitudinal (ny) 
and transverse {v±) velocities are calculated using a lattice 
spacing of 3.8A and hopping t = OAeV. 



transition between hole- like and electron- like. 

The sign of the Hall coefficient depends sensitively on 
the curvature of Fermi surface and the Fermi velocity 
and small changes can switch the sign. The fact that 
experimentally the sign of the Hall coefficient in strong 
magnetic fields appears to be non-universal may be an 
indication that it is related to such open orbits. A more 
realistic calculation of the stripe band-structure is needed 
including also an examination of the total Hall coefficient 
due to a possibly electron-like open Fermi surface section 
together with hole pockets. 
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FIG. 3: (Color online) Spectral intensity at the Fermi energy 

(a) together with the local density of states of the antinodal 

(b) and nodal (c) states plotted on the same normalized scale 
with respect to the max;imum density (black). Note that the 
high density regions are shifted by two sites between nodal 
and antinodal weight. The antinodal weight has high density 
at the AF domain walls of the stripe potential whereas the 
nodal weight has high intensity between. Also shown is the 
amplitude of the Fourier transform of the LDOS transverse 
to the stripe direction with the uniform density component 
k± = normalized to one. The LDOS modulations in both 
regions is roughly 30% of the mean density. (The combined 
total LDOS modulation at the Fermi energy, not displayed, is 
only around 5%. The integrated charge modulation over all 
occupied states is around 4%.) 



Charge order 

Let us also consider the local density of states of this 
system LDOS{x,cj) = \^^{x)\^5{uj - tu^). We will 
be interested in considering the connection between the 
spectral function and the local density of states at the 
Fermi energy and in particular the distinction between 
nodal and antinodal states. Since momentum is con- 
served modulo TT along the stripe direction fcy we can 
associate to every state a momentum — 7r/2 < k\\^a < 7r/2 
and choose to make a partial sum in the LDOS that only 
includes states with certain fcy . We define 

LDOSnodal = \^a{x)\'^5{uJ-UJa) 
Q:|fe||,„|>Tr/4 

LDOSa-nodal = \^ c.{x)\H{u - U^) (2) 

Q:|fc|l j,|<ir/4 

Figure [3] shows the nodal and antinodal LDOS for the 
period four array, the antinodal states have a higher den- 
sity on sites x mod 4 =0 or 1 while the nodal states have 
a higher density on sites x mod 4 =2 or 3 (with respect to 



the potential V4^o(a^ mod 8) = 1, 1, —1, 1, —1, —1, 1, -1). 
As emphasized in several other publication^2llll] the 
nodal states derive from the antiferromagnet whereas the 
antinodal states are the "stripe states" living on domain 
walls of the antiferromagnet. We will find that this dis- 
tinction will be crucial when we in the next section con- 
sider a system without long range spin stripe order. 

IV. DISORDERED STRIPES 

We will now consider a charge ordered but spin disor- 
dered system. We generate a stochastic spin-stripe po- 
tential with a short correlation length through an up- 
date procedure V{x + 1) = Vq{x)Vo{x + l)V{x) + 
where —2m < i/i < 2m is a random number and ap- 
ply a cutoff —m < V(x) < m at every step. Here 
Vo{x) gives the sign of the potential for the correspond- 
ing spin-stripe ordered array which in the case of the 
simple charge-period four bond-centered stripe takes the 
value Vo{x mod 8) = 1, 1,-1, 1, -1,-1, 1, -1. The dis- 
tribution is readily evaluated analytically and we find a 








FIG. 4: (Color online) Disorder averaged spin density corre- 
lator transverse to the stripe extension for period four bond- 
centered stripes showing the approximate self-consistency of 
the stripe potential. In (a) the absolute value in real space and 
in (b) the Fourier transform. Solid points are values for the 
stripe potential and open symbols for the actual spin densities 
which follow from the potential. The latter are multiplied by 
= (3.5)^ as discussed in the text. In (a) the position is 
on a stripe (triangles) or between stripes (stars), where the 
smaller amplitude on a stripe is consistent with a higher hole 
density. The lines are analytic results for the disorder aver- 
aged potential, Equations |3] (absolute value and m = 1) and 
|4]in (a) and (b) respectively. The average is taken over 5000 
runs with system size 120x120 (unit cell 120x2). 



disorder averaged correlation function 

{Vix)V{x + j)) = '^m^Vo{x)V„ix + j)e-''^^^ , (3) 

with a short correlation length ^ = 1/k — 1/ ln(2) w 1 .44. 
We will take m = 1 which we find gives a good correspon- 
dence with the ordered system considered previously, i.e. 
the same Fermi energy gives roughly the same doping. 
In order to be able to consider large system sizes (in the 
transverse direction) we still assume antiferromagnetic 
order along the stripes. Assuming short-range correla- 
tions would broaden the states in fcy but as far as the 
transverse charge order is concerned we do not expect 
this to have any significant impact. For the numerical 
calculations we use a transverse system size N which is 
large compared to the correlation length of the poten- 
tial (typically N — 120) and diagonalize the system for 
a given static but random configuration of the potential 
and take the average over many realizations. 

Physically what we have in mind is a system with fiuc- 
tuating short-range spin stripe correlations with some 
correlation length ^ and frequency luq . The disorder aver- 
age can thus be relevant either for a fast, compared to ujq, 
but extended probe such as inelastic neutron scattering 
or ARPES (disorder average in real space) or a slow but 
local probe such as STS (disorder average in time). There 
is of course a significant approximation involved in repre- 
senting the spin fluctuations by a sequence of static spin 
configurations. In particular quasi-particle dynamics will 
not be accurately represented using such a quasi-static 
spin approximation. However for the gross and qualita- 
tive features of the spectral weight distribution we expect 
that the approximation can give reasonable results. 
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FIG. 5: (Color online) Spectral intensity (a) of spin disor- 
dered system (same as in Fig. |4| together with the LDOS at 
the Fermi energy of antinodal (b) and nodal states (c). The 
LDOS is shown for an arbitrary section of the system and is 
plotted on the same relative scale as for the spin-stripe or- 
dered system (Fig. [3|. The relative variations is around 18% 
for the antinodal LDOS and around 3% for the nodal LDOS. 
The components at momenta n/2 and 37r/2 correspond to 
roughly 15% density variations of the antinodal LDOS (d) 
and 1.5% (barely visible in e) for the nodal LDOS variations 
with the rest of the variations being due to random disorder. 
The total density modulation integrated over all filled states 
is around 3%. The simulation is averaged over 10000 runs 
with system size 120x120. 



Spin correlations 

We think of the stripe potential as reflecting the in- 
stantaneous spin configuration. Clearly this disordered 
potential will not be a self-consistent mean field solu- 
tion to the Hubbard model. Nevertheless, to be phys- 
ically reasonable we may at least require that the dis- 
order averaged potential should be roughly self consis- 
tent in the sense that it faithfully reflects the average 
spin conflguration. For the potential as we have cho- 
sen it every site is equivalent up to signs given by Vq, 
however for the actual disorder averaged spin density 
we need to distinguish between two inequivalent sites 
which are on average on a stripe {x mod 4 =0 or 1) 
or between a stripe {x mod 4 =2 or 3) with respect 
to the potential Vo{x mod 8) = 1, 1, —I, I, —I, —1, 1, — 1. 
Figure [4] shows the corresponding spin-spin correlators 
\{S^{x)S^{x + as well as the Fourier transforms 

Spotiq) = {{l/N)j:^^^^^V{xi)V{x2)e^'i(^'--''>) and 

S{q) = ((l/iV)^^^_^^5^(a;i)^^(a;2)e*'?("^-"i)). From 
the analytic expression for the correlation function of the 
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potential, Eqn. [3] with m = 1, we can calculate 

2/ 2^x/2 e'' - cos(g ~ gp) 

i>pot(q) - 3^ ^+ 8 2^ cosh(/«) - cos(g - go) 

qo=±-l^/i 

^ 2 + 72 e'^ - cos(g ~ gp) \ 

8 ^ COsh(K) - COs(g - On)'' 

with K = ln(2). The spin-spin correlators are multi- 
plied by a factor — (3.5)^ and the good agreement 
with the potential indicates that the averaged potential 
is quite close to being self-consistent for the Hubbard 
model with U = 3.5. The main source of the discrep- 
ancy is caused by the fact that the potential assumes 
on average a uniform magnitude of the local magnetiza- 
tion whereas a self-consistent calculation gives a smaller 
magnetization on the stripes (Fig. |4^) consistent with a 
larger hole density. This can be fine tuned without any 
significant impact on the subsequent results and we will 
prefer to keep the model as simple as possible. 



Charge order 

Let us now study the charge density in this system 
with short range spin correlations. Here we consider 
two objects, the disorder averaged local density of states 
LDOS{x,uj) = {LDOS{x,uj)) and the disorder averaged 
spectral function A(k, lo) — {A{k, ut)) where we can again 
divide the LDOS into nodal and antinodal parts accord- 
ing to Eqn. [2] 

Figure [5] shows the spectral weight integrated over a 
window {5w = 0.05)ES1 close to the Fermi energy, as well 
as the corresponding nodal and anti- nodal LDOS. We 
take the Fermi energy the same {Ep = —1.15) as for the 
period four stripe ordered system considered above which 
corresponds here to an average hole doping of around 
18%. Due to the disorder the spectral weight is spread 
out such that the stripe Fermi surface and nodal pocket 
loose their integrity and merge into a single large Fermi 
surface which if we symmetrize with respect to stripes 
along both directions would close around the (tt, tt) point. 

What is particularly notable are the results for the 
LDOS where the nodal LDOS is almost uniform with 
only a very weak signature of the charge order whereas 
the antinodal LDOS shows long range charge-stripe order 
not very different in magnitude from that of the stripe or- 
dered system considered previously. For the nodal states 
there is an apparent restoration of the rotational sym- 
metry of the lattice for the spin disordered system. The 
results can be roughly understood based on the knowl- 
edge of the stripe ordered system in which the antinodal 
weight comes from stripe states that are centered on an- 
tiferromagnetic domain walls and are quite insensitive to 
whether or not there is long-range and static spin order. 
The nodal states on the other hand live primarily be- 
tween stripes for a spin-stripe ordered state and the spin 



fluctuations apparently conceal the stripes. A comple- 
mentary understanding can be found from earlier work 
on the same model where it was noted that in the pres- 
ence of weak potential disorder the localization length 
transverse to the stripes of the antinodal states is of the 
order of or less than the stripe period whereas the local- 
ization length of the nodal states can be several times the 
stripe period.-^ 



V. DISCUSSION 

The Fermi momentum resolved charge order found 
here is yet another manifestation of what has been called 
the "nodal-antinodal dichotomy" in which the low-energy 
spectral properties are very different for the nodal and 
antinodal regions.— In particular the enigmatic pseudo- 
gap is known from ARPES to be exactly such a mo- 
mentum resolved phenomenon in which there is a partial 
suppression of the antinodal spectral weight. 

Let us explore the implications of an antinodal pseu- 
dogap on our results in a qualitative fashion: First, it 
would lift the bulk of the antinodal spectral weight away 
from the Fermi energy leaving only the nodal Fermi arc. 
In addition it would imply that the charge-stripe order 
may only be visible above the pseudogap energy scale. 
The latter property does indeed seem to be realized for 
the "stripe glass" observed by tunneling spectroscopy for 
which the differential tunneling conductance is nearly 
uniform at low energies and only show the unidirecti onal 
modulations at energies around the pseudogap energy.'^'^ 
These observations together with the results presented 
here give a very strong indication that the pseudogap is 
in fact related to stripes. Also corroborating this picture 
is ARPES data which find that the pseudogap disappears 
abruptly with the momentum along th e Fermi surface in 
going from (7r,0) to (7r/2, 7r/2^p2lilEil in a fashion analo- 
gous to the charge stripe order presented here. 

The pseudogap is then a strong correlation gap (spin- 
gap) on a single stripe characteri zed as an interacting 
quasi-one dimensional electron gas.^^^ In the work pre- 
sented here any strong correlation effects apart from the 
stripe order itself are completely ignored but we do find 
that the antinodal spectral weight comes from the quasi- 
one dimensional stripe states thus providing a logical con- 
sistency with an antinodal pseudogap from stripes. One 
way to argue why it may be reasonable to treat the stripe 
correlations separately is if there is a separation of energy 
scales in the renormalization group sense in which the 
antiferromagnetic correlations, which for a doped system 
correspond to stripe correlations, dominate and residual 
correlations which gives rise to t he pse udogap become 
prominent at a lower energy scale .'^'^^EIl 

As is well known for the underdoped cuprates the 
superconducting order occurs at an even lower energy 
scale. F rom, for example, ARPES results on the gap 
evolutioiJi^Ill! and STS results on the Bogoliubov quasi- 
particle interference it is by now quite clear that long- 
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range superconducting order is intimately linked to the 
nodal liquid which as we have found can be isotropic even 
in the presence of charge stripe order. In this context 
we note that the spin stripe ordered Lai.875Bao.i25Cu04 
has a significantly suppressed superconducting transi- 
tion temp eratur e while there is still a large antinodal 
pseudogap ! ^ ^ 1"^^ I The suppression of superconductivity in 
this system may thus be related to anisotropy of the 
nodal Uquid due to spin stripe order 

As a continuation of this work the modeling of dis- 
ordered stripes should be extended to include also some 
characterization of the pseudogap and a d-wave like pair- 
field which could allow for a more detailed comparison 
with STS and ARPES results. For the ordered stripe 
phases work is in progress to calculate the total conduc- 
tivity resulting from the contributions of hole pocket and 
open stripe band and to characterize more generally the 
sign and oscillations of the Hall resistivity in such a sys- 
tem. 

To summarize, we find that for a charge ordered but 



spin disordered stripe state there is a natural separation 
into nodal and antinodal states for which the latter shows 
a much stronger signature of the charge order. The local 
density of states at the Fermi energy is almost homoge- 
neous when integrated only over states that contribute to 
the nodal region spectral weight. We also consider com- 
mensurate spin and charge stripe ordered systems where 
the Fermi surface consists of a nodal hole pocket and 
an open (quasi-one dimensional) antinodal section. The 
relation between the size of the pocket and the hole den- 
sity depends on whether the charge periodicity is odd or 
even and we give upper estimates of the possible size of 
a pocket for charge period four and five of 1.6% and 4% 
respectively. We find that the open section may be elec- 
tron or hole like depending on the details of the stripe 
band structure. 

We thank John Tranquada and Steven Kivelson for 
comments and Seamus Davis for discussions of his exper- 
imental dataP^. 
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